Persistent currents and quantised vortices in a polariton superfiuid 
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Semiconductor microcavity polaritons in the 
optical parametric scattering regime have been 
recently demonstrated to display a new variety 
of dissipationless superfiuid behaviour. We re- 
port the first observation in resonantly pumped 
exciton polaritons of a metastable persistent su- 
perflow carrying quantum of angular momentum, 
to. The quantised vortex, excited by a weak 2ps 
pulsed probe, is shown to last for at least 80ps, 
limited only by the leaking outside the cavity. 
The polariton circulating superfiuid persists in 
the absence of the driving rotating probe with 
no apparent dissipation. In addition, for a mov- 
ing superfiuid, we show the coherent splitting of 
a quantised double vortex, with charge to = 2, 
into two singly quantised vortices of m = 1. Re- 
markably, we observe the to = 2 vortex to be sta- 
ble when they are at rest. The experimental re- 
sults are compared with a theoretical analysis, ob- 
tained describing the triggered parametric scat- 
tering regime of polaritons via a two-component 
Gross-Pitaevskii equation, including pump and 
decay processes. 

In the past few decades there has been a strenuous 
search for macroscopic coherence and phenomena related 
to Bose-Einstein condensation (BEC) in the solid state. 
The first realisation of a BEC in semiconductor microcav- 
ities [H has inaugurated a new era in the study of strongly 
coupled light-matter systems. The growing interest in 
this field can be attributed to the unique properties of 
exciton polaritons in microcavities the composite 
particles resulting from strong light-matter interactions. 
Properties of a polariton condensate Q differ from that of 
other known condensates, such as ultracold atomic BECs 
and superfiuid 4 Hc. In particular, polaritons have a short 
life time of the order of picoseconds, therefore needing 
continuous pumping to balance decay and reach a steady 
state regime. Rather than a drawback, the intrinsic non- 
equilibrium nature enriches the features of polariton con- 
densation, but at the same time poses fundamental ques- 
tions about the robustness of the coherence phenomena 



to dissipation and non-equilibrium. Superfiuid proper- 
ties of non-equilibrium condensates in dissipative envi- 
ronment still need to be understood [i| and this letter 
provides a substantial advance in this direction. 

One route to experimentally inject polaritons into a 
microcavity is by non-resonant (incoherent) pumping. 
For incoherent pumping, polaritons have been shown 
to enter, within their short lifetime, a macroscopically 
coherent BEC phase [H, [f|. As a consequence of non- 
equilibrium currents in an inhomogeneous system, carry- 
ing particles from gain to loss dominated regions, deter- 
ministic quantised vortices have been observed to spon- 
taneously appear Q. However, the unusual — diffusive at 
small momenta — form of the excitation spectrum 0, S] , 
which hinders the fulfilling of the Landau criterion, puts 
under debate the possibility of dissipationless superflow 
in incoherently pumped polariton systems. 

A different scenario characterises coherent resonant in- 
jection of parametrically pumped polariton condensates, 
which have been recently shown to exhibit a new form 
of non-equilibrium superfiuid behaviour [Tell . In the 
optical parametric oscillator (OPO) regime [ll|, bosonic 
final state stimulation causes polariton pairs to coher- 
ently scatter from the pump state to the signal and idler 
states, which, at threshold, have a state occupancy of or- 
der one. The properties of the quantum fluids generated 
by OPO at idler and signal have been recently tested via 
a triggered optical parametric oscillator (TOPO) config- 
uration Q . An additional weak pulsed probe laser beam 
has been used to create a traveling, long-living, coher- 
ent polaritons signal, continuously fed by the OPO. The 
traveling signal has been shown to display superfiuid be- 
haviour through linear dispersion and frictionless flow. 
However, generating long-lived quantised vortices and 
the possibility of metastable persistent flow in this con- 
figuration still remain missing in the superfiuid 'check- 
list' 

We show for the first time that the non-equilibrium su- 
perfiuid generated with a TOPO can hold a vortex state 
carrying angular momentum with quantum number to, 
and this is observed up to 80ps after the probe has been 
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switched off, more than 10 times the lifetime of photons 
in the microcavity. The vortex is excited by a weak, 
pulsed probe resonant with the signal (see Fig. []} lasting 
2ps only. In other words we stir the polariton superfluid 
only for a short time and observe its long lived rotation 
with the identical quantum of angular momentum on a 
time scale almost 40 times longer than the duration of 
the pulse (see Figs. [2] and [3]), the analogous of the rotat- 
ing drive in trapped gases. Thus in this work, we observe 
two of the most fundamental properties of a superfluid: 
Quantisation of circulation and stability of persistent cur- 
rents in the novel non-equilibrium and dissipative system 
of driven polariton condensates, demonstrating robust- 
ness of superfluid phenomena. 

Another fundamental property of vortices with higher 
charge number m observed in atomic superfluids, is the 
tendency to coherently split into many vortices of charge 
m = 1. To study the stability properties of vortices in 
polariton condensates, we have injected a vortex with 
m = 2 and let it evolve within the time of observation 
(see Fig. 2]). Surprisingly, we found that the vortex splits 
only when it is moving. Indeed, dynamic instability of 
doubly quantised vortices, in atomic BEC, has been ob- 
served [12J, where a double quantised vortex has been 
showed to coherently break into two singly quantised vor- 
tices. However, there are few systems where double quan- 
tised vortices have been shown to be stable, for example 
in superconductors in the presence of pinning forces [1 31 ] , 
and in superfluid 3 He-A, which has a multicomponent 
order parameter 
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Remarkable stability of m = 2 
static vortex in polariton systems provides an additional 
example. 

In our experiment we use a semiconductor microcavity 
with a Rabi splitting = 4.4meV and the cavity pho- 
ton energy slightly negatively detuned below the exci- 
ton energy, S E [—3, — l]meV — for details see the Method 
section — . A Ti-Sapphire laser is tuned in resonance with 
the lower polariton branch (LPB), injecting polaritons 
close to the point of inflection and giving rise to a contin- 
uously pumped OPO (sec Fig. [TJ) . Above a pump thresh- 
old, the generated signal state close to zero momentum 
(k s = 0) as well as the idler state at high momentum, 
form an out of equilibrium coherent polariton superfluid. 
At a given time we trigger a new scattering process on 
top of the OPO (TOPO) signal with a resonant pulsed 
probe smaller than the signal. The probe is a pulsed 
Laguerre-Gauss beam carrying a vortex of given angu- 
lar momentum m, the phase of which winds around the 
vortex core with values from to 2-nm. However, after 
a short time, ~ 2ps, we remove the probe, leaving the 
polariton coherent state free to rotate, without the driv- 
ing field. While even a classical fluid acquires angular 
momentum in the presence of an external rotating drive, 
only a superfluid can exhibit long lived circulating flow 
once the external drive is turned off. 

In order to demonstrate persistence of the vortex an- 
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FIG. 1: Lower (LPB) and upper polariton branch (UPB) dis- 
persions together with the schematic representation of the 
TOPO excitation (a). Resonantly pumping the LPB initi- 
ates, above threshold, stimulated scattering to a signal close 
to zero momentum and an idler at higher momentum. The 
m = 1 Laguerre-Gauss beam resonant with the signal (b) is 
used as a weak pulsed triggering probe to stir the superfluid. 
The corresponding interference image is shown in panel (c). 



gular momentum, we detect the phase pattern generated 
by making interfere the TOPO signal with an expanded 
and flipped spatial region far from the vortex core (where 
the phase is approximately constant) in a Michelson in- 
terferometer. A fork-like dislocation with a difference of 
m arms corresponds to phase winding by 2iTm around 
the vortex core (see lower panel of Fig. [I}. 

Using a streak camera we can follow the evolution in 
time of the vortex generated by the pulsed probe. Time 
intervals of 1 or 2ps are used to reconstruct 2D images of 
the signal state after the perturbation has arrived. Every 
picture is the result of an average over many shots, all 
taken at the same time and same conditions. In order 
to separate the contribution of the signal from that of 
the pump, both in theory and experiments, we filter the 
signal images in momentum space in a cone around k = 
of approximately ±7°. 

The time evolution of an excited m = 1 vortex to- 
gether with its interference pattern, which characterises 
unequivocally the vortex state, are shown in Figure [5] 
The probe has a power of 0.3mW which is comparable 
to the OPO signal emission intensity, before the probe 
arrives on the sample at time t p b- The effect of the probe 
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FIG. 2: Time evolution of the TOPO signal after a weak 
pulsed probe with a vortex of m = 1 has been excited (d- 
f). The corresponding interference images (a-c) are obtained 
by overlapping the vortex with a small expanded region of 
the same image far from the vortex core, where the phase 
is constant. The time origin is taken when the TOPO has 
reached 80% of its maximum value. 
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FIG. 3: Numerical simulations of the time evolution of the 
TOPO signal after a weak pulsed probe with a vortex of m — 1 
has been impinged at k = 0, and for pumping strength very 
close to the OPO threshold (/„ = 0.9994 th) ). (a) and (b) 
show images of the signal obtained by momentum filtering in 
a cone of around 7°, while interference fringes of (c), showing 
the characteristic fork-like dislocation of the m = 1 vortex, 
are obtained by taking the full emission without filtering out 
the pump. 



is to enhance the polariton signal emission by a factor 
of 20 with a delay of lOps after the pulse arrival time. 
The vortex appears as a new state on top of the OPO 
signal. To eliminate the contribution of the steady state, 
we subtract the emission of the OPO in absence of the 
probe pulse. The signal generated by this perturbation 
lasts about a hundred picoseconds. 

The images of Fig. [5] demonstrate that the effect of 
the probe on the signal steady state is lasting for at least 
80 picoseconds after the probe is switched off, and during 
this time we observe a circular flow of polaritons rotating 
unperturbed. This effect is visible in the real space im- 
ages through the 'thoroidal shape' of the polariton emis- 
sion which is present at any time. During the time the 
vortex lasts, we can observe the vortex slowly drifting to 



the right and changing slightly its shape. The drift is 
due to the fact that the probe has a small component 
at finite positive momentum. The small changes in time 
of the polariton vortex shape are due to the presence 
of local defects 
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which tends to influence the polari- 
ton density distribution around the vortex core, affecting 
therefore the shape of the excited vortex. However, the 
interference pattern clearly shows that the spatial phase 
relation, including the fork dislocation, remains constant 
and with full contrast all over a wide area, even for times 
when the interference contrast of the signal polariton has 
dropped down to 30% due to the relatively short (50ps) 
polariton coherence time — see figure in the supplemen- 
tary material — . This means that even when polaritons 
have lost any phase relation in time, the angular momen- 
tum is still conserved. 

The sequence in Figure [2] demonstrates that the vor- 
tex remains steady as a persisting metastable state as 
long as the extra population, created by the probe, lasts. 
Eventually the system will tend to the steady state and 
all angular momentum will be transferred to the photons 
leaking out of the cavity. However, down to the last de- 
tectable signal we observe a decrease of vortex radius and 
intensity but no loss of winding number due to dissipa- 
tion. This is revealed by the strong contrast of the fork in 
the interference images. In the opposite scenario, the in- 
terference images would show low contrast in the vortex 
core region, indicating a mixture of the population which 
has undergone dissipation and the one still carrying the 
quantised angular momentum. 

We have numerically simulated the experiments 
by making use of mean-field two-component Gross- 
Pitaevskii equations for the coupled cavity and exciton 
fields with external pump and decay — see also Method 
section — . In order to reproduce conditions similar to the 
experimental ones (TOPO), where the vortex appears as 
a new state on top of the OPO signal, we choose parame- 
ters very close to the OPO threshold from below. Under 
these conditions we observe the weak pulsed probe to 
trigger a vortex state at the signal which can last around 
80ps (sec Fig. [3]) , in agreement with the experimental re- 
sults. Images of the vortex are obtained by filtering in 
momentum space. Interference fringes showing the char- 
acteristic fork-like dislocation of the m = 1 vortex are 
obtained by taking the full emission without filtering out 
the pump. Like in the experiments, we observe the vor- 
tex state lasting as long as the signal lasts. In contrast, 
a vortex created at k = well below the OPO thresh- 
old, typically decays already two orders of magnitude in 
intensity within the firsts 20ps. 

A second experiment aimed at investigating the stabil- 
ity of doubly charged vortices in different regimes. There- 
fore, we repeat both our experiments and simulations in 
conditions similar to the ones previously described, but 
now generating a TOPO doubly quantised vortex, i.e. 
with m = 2. By repeating the experiment in a variety 
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FIG. 4: Time evolution of a m = 2 vortex TOPO signal ex- 
cited close to zero momentum, |k| = 0, (a-c) for which neither 
motion nor splitting of the vortex could be detected. However 
for a signal at k s pointing to the right (d-f), we can detect 
the m — 2 vortex both moving and splitting into two single 
quantised m = 1 vortices. Numerical simulations of the time 
evolution of the TOPO signal in the same conditions of Fig. [3] 
(g-i).The doubly quantised vortex splits already after 6ps the 
probe has been switched off, corresponding to the moment at 
which the vortex starts also drifting. In the last panel of each 
row we plot the interference images corresponding to a late 
time. 



of different conditions, we observe that when the vortex 
is excited on a static signal centered at k = the vor- 
tex does not split within its lifetime. However, exciting 
the TOPO signal with a finite momentum, thus making 
it moving inside the pump spot, we observe the doubly 
quantised vortex splitting into two singly quantised vor- 
tices as shown in Fig. 2] Our theoretical analysis gives 
the same result: vortex at rest is stable while a moving 
vortex splits (see lower panels in Fig. [4]) , obtained under 
the same conditions as in Fig. [3J Splitting occurs when 
the vortex starts slowly drifting in space and the signal 
emission becomes asymmetric around k = 0, which hap- 
pens around 6ps after the probe has been switched off. 
The splitting of the vortex is confirmed by the interfer- 
ence image. 

In the numerical simulations, we also observed a quali- 
tatively different regime in which very long lived vortices 
are found to appear in the steady state signal, following 
the probe pulse. In the simulation of Fig. once the 



OPO steady state is reached, we turn on the probe reso- 
nant with the OPO signal for 2ps only. After the probe 
is switched off, the signal experiences a transient period, 
when either the vortex drifts around or no vortex can be 
detected. After the transient, independently on the in- 
tensity and size of the probe, a vortex with a quantum of 
angular momentum stabilises into the signal and lasts as 
long as our simulations (~ Ins). We have observed sta- 
ble vortex solutions for probe intensities down to 1/20 
the intensity of the signal and also for different vortex 
probe sizes. A stable vortex solution does not always ex- 
ist; there are cases where, during the transient period, 
the excited vortex spirals out of the signal. Stable vortex 
solutions above the OPO threshold were also reported in 
Ref. [Tj|. 

By contrast, in experiments, vortices in the steady 
state regime are not detected. One possibility could be 
that disorder can pin the phase of the signal OPO, thus 
making it harder for the probe to change the signal phase. 
A second reason could be that the presence of random 
noise (temperature and laser power fluctuations, normal 
fluid component, etc.) can contribute to give different 
dynamical scenarios at any experimental shot, depend- 
ing on the initial conditions. Indeed, differently from the 
simulations, the experimental images are obtained by ac- 
cumulating over several million realizations of the exper- 
iment, hiding the different final spatial configuration of 
the stable vortex in the signal OPO. However, in case of 
the TOPO configuration, the experiments start always 
with the same initial condition, which is imposed by the 
probe and kept by the driving field, and there is no tran- 
sient period. Moreover, the presence of a strong OPO 
signal, which fills the sample inhomogeneities, strongly 
contributes in smoothing local defects, helping the free 
rotation of the triggered polariton condensate. 

To conclude, we have demonstrated the presence of 
long lived quantised currents in a polariton supcrfluid 
which last as long as polaritons are present in that state. 
Our report is the first observation of persistent currents 
in semiconductors. The superfluid phenomena observed 
in our microcavity is shown to be robust despite the rel- 
atively large sample temperature (10K) and the presence 
of strong decay and dissipation. Moreover, the unex- 
pected stability of winding number m = 2 static vortices 
observed in this work, as well as our theoretical predic- 
tion of having, above OPO threshold, an imprinting into 
the OPO signal of the winding number of the probe, in- 
dicate the richness of superfluid behaviours which can be 
realised with driven polariton condensates. 



METHODS SUMMARY 

In our TOPO experiment Q , the pump power was kept 
always above threshold for the OPO, in a pumping regime 
where the energy rcnormalisation around the spot area 
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FIG. 5: Simulations showing the persistent vortex supercur- 
rent with one quantum of angular momentum, m = 1, im- 
posed with a short, pulsed vortex probe at the OPO signal, 
for pump strength above threshold f p — 1.12/p . Simula- 
tions show the existence of a stable vortex solution even ~ Ins 
after the driving probe has been removed. 



due to interaction is as constant as possible. The vortex 
probe state is prepared by scattering a pulsed laser Gaus- 
sian beam in an hologram with a single or double fork- 
like dislocation. This gives rise to a first order Laguerre- 
Gauss beam with a winding number respectively cither 
to = 1 or to = 2. The probe beam is focused at the center 
of the cw pump in resonance with the signal emission en- 
ergy and at around k = 0. Every picture is the result of 
an average over many shots, as single shot measurements 
would give a too low signal to noise ratio. 

The dynamics of amplitudes and phases of TOPO is 
computed with a two-component Gross-Pitaevskii equa- 
tion with external pumping and decay for the coupled 
cavity and exciton fields ipc,x( r ,t) (fi = 1) : 
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Since the exciton mass is four orders of magnitude larger 
than the photon mass, we neglect the excitonic dispersion 
and assume a quadratic dispersion for the cavity photon, 
loc = ^c(0) — 2^:- The fields decay with rates nc.x- 
f2# is the photon-exciton coupling. The cavity field is 
driven by an external cw pump field, 
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while the probe is a Laguerre-Gaussian pulsed beam, 
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producing a vortex at r p b with winding number to . The 
exciton repulsive interaction strength gx can be set to 
one by rescaling fields and pump strengths. We solve 
Eq. |T]) numerically by using the 5 th -ordcr adaptive-step 
Rungc-Kutta algorithm on a 2D grid. 



METHODS 
Experiments 

The sample studied is a A/2 AlAs microcavity with 
a 20 nm GaAs quantum well placed at the antinode of 
the cavity electromagnetic field. The cavity is formed 
by two high reflective Bragg mirrors of 25 pairs at the 
bottom and 15.5 on the top of the structure. The pump is 
obtained with a cw Ti-Sapphire laser, resonantly exciting 
the LPB close to the inflection point at 9° . The spot size 
is of ~ 150/im. A pulsed Ti-Sapphire laser excites a probe 
2ps long focussed to a diameter sa 250/im. 

The experiments are performed at a cryogenic temper- 
ature of 10 K and using a very high numerical aperture 
lens (0.6) so that the sample could be accessed by angles 
as large as 25° and the photolumincsccncc (PL) can be si- 
multaneously collected from the signal state in the near 
as well as the far field. The PL was collected through 
a 0.5m spectrometer into a streak camera, working in 
synchroscan mode, with 4ps time resolution, allowing for 
energy- as well as time-resolved images. In order to reach 
high time resolution, all the images where obtained by 
filtering the signal in the far field without spectrally re- 
solving the emission of the OPO states. 



Theory 

In our numerical simulations, we take mc = 2.3 x 
10 _5 too and a Rabi splitting — 4.4meV determined 
experimentally, fixing the detuning to zero, S = wc(0) — 
(0) = 0. In the regime considered by the experiments, 
we can neglect the saturation of the dipole coupling [ItJ • 
We pump at ~ 12° (corresponding to |k p | = 1.6/im _1 for 
a LPB energy at normal incidence of wlp(O) = 1.524eV). 
The energy of the pump io p is tuned around 0.5meV 
above the bare LPB, in order to take into account the 
energy renormalisation due to interaction We have 
chosen a pump profile with full width at half maximum 
(FWHM) of 35/im. 

The pump strength f p can be chosen to be below or 
above threshold for OPO. In particular, we consider two 
different situations in our numerical simulations. In one 
case we fix the pump strength f p very close to the thresh- 
old for OPO, f p = 0.999/i th) . We switch on for 2ps only 
a probe with a to = 1, 2 vortex with a p \, = 4.5/im, mo- 
mentum k p b = and energy u) p b around 0.5meV above 
the bare LPB dispersion. The probe triggers a signal, 
which intensity first quickly grows and then slowly de- 
cays in time — after 60ps the signal intensity is reduced 
by 2 orders of magnitude — . The vortex excited by the 
probe starts very slowly drifting inside the signal after 
6ps and is visible for around 80ps. 

In the second class of numerical simulations, we drive 
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the system into the OPO regime at f p = l.l2fp th \ 
which generates a signal narrower than the pump with 
FWHM=12/im at slightly negative angles (~ -2°). The 
signal is shifted parallel to the direction of the incident 
pumping laser of about 5^tm [l5| ■ Once the steady state 
is reached for the pump only OPO configuration, we turn 
on the probe (Equation ([2|)) for 2ps only. The momentum 
k p b of the probe is tuned to the one of the OPO signal 
and its energy set around 0.5meV above the bare disper- 
sion. We choose a probe with a vortex of m = 1. We 
have considered different sizes of the probe vortex, from 
2c p h = 9^tm down to 1/xm. In addition, we shift the vor- 
tex core center to r p b = (4, 0)/im, in order to overlap the 
core with the maximum signal intensity. For the probe 
intensity we have considered different values in the range 
[1/200, 1] times the intensity of the pump, corresponding 
to a range of [1/20, 10] times the intensity of the signal. 
We have checked that neither the intensity nor the size 
of the probe changes qualitatively the outcome of the 
simulations. The existence of a stable vortex solution is 
independent on both the probe intensity and the probe 
size. Remarkably the signal rearranges itself during the 
transient period and generate a stable vortex configura- 
tion with a vortex core size independent on the one of the 
probe. The transient period can last between [50, 300]ps, 
depending on the probe intensity and size. Our conclu- 
sion is that the probe can only change the quantitative 
behaviour and length of the transient period. 
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The calculated images of the signal in Figs. Efa),(b), 
|4jg),(h) and[Ua),(b) are obtained by filtering in momen- 
tum space around the signal momentum, k s ~ 0. We set 
to zero all the emission in momentum space aside the one 
coming from the signal and fast Fourier transform back 
to real space. In this way the strong emission coming 
at the pump angles is masked out. Interference fringes 
instead are obtained by taking the full emission without 
filtering out the pump (and the idler). In numerical sim- 
ulations we consider a single shot measurement at a given 
time, therefore pump, signal and idler give rise to inter- 
ference fringes, even if they oscillate at different energies. 
Such interferences are clearly washed out in experiments, 
where single shot measurements are hindered by the low 
signal to noise ratio. For this reason, in experiments we 
need to use a Michclson interferometer. 
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SUPPLEMENTARY MATERIAL 



In Fig. [6] we show the intensity of the interference 
fringes of the vortex with m = 1 (Fig. [2] of the main 
text) obtained with the Michelson interferometer but at 
different time delays between the two arms. Immedi- 
ately after the probe has arrived, and at zero delay, the 
fringes exhibits a visibility close to 100%, while it reduces 
to 30% after a delay of 55ps. However the interferences 
of the vortex core of Fig. [2] — taken always at zero de- 
lay time — do not show any degradation of the visibility 
even 80ps after the probe has arrived, demonstrating the 
persistence of the polariton vorticity. 
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FIG. 6: Intensity of interference fringes (arbitrary units) of 
the polariton emission after arrival of the probe pulse. The 
two images are taken at different delay times between the 
polariton signal and its mirrored counterpart: (a) Ops and 
(b) 55ps time delay. 



